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Abstract: We suppose that the baryon asymmetry is produced by thermal leptogenesis 
(with flavour effects), at temperatures ~ 10 9 — 10 10 GeV, in the supersymmetric seesaw with 
universal and real soft terms. The parameter space is restricted by assuming that ► 
processes will be seen in upcoming experiments. We study the sensitivity of the baryon 
asymmetry to the phases of the lepton mixing matrix, and find that leptogenesis can work 
for any value of the phases. We also estimate the contribution to the electric dipole moment 
of the electron, arising from the seesaw, and find that it is (just) beyond the sensitivity of 
next generation experiments 10 _29 e cm). The fourteen dimensional parameter space is 
efficiently explored with a Monte Carlo Markov Chain, which concentrates on the regions 
of interest. 
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1. Introduction 

Neutrino masses are evidence for beyond the Standard Model (SM) physics. A simple 
extension of the standard model that accounts for neutrino masses is the seesaw mechanism 
[1], where heavy majorana right-handed neutrinos are added to the SM. Moreover, the 
seesaw scenario provides a very attractive framework to explain the baryon asymmetry 
of the universe (BAU) through the leptogenesis [2] mechanism, without inducing proton 
decay. 

CP violation is a necessary ingredient to explain the BAU and, if this asymmetry 
is produced via leptogenesis, the required CP violation is encoded in the CP violating 
phases of the lepton sector. Three of them are the well known Dirac and Majorana phases 
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of the PMNS mixing matrix, that are in principle measurable. Any observation of CP 
violation in the lepton sector, for instance CP violation in neutrino oscillations due to 
the PMNS phase <5, would then support leptogenesis by demonstrating that CP is not a 
symmetry of leptons. However, even in this very promising case, the question of whether 
the BAU is produced via leptogenesis is far from being answered, because it is not possible 
to reconstruct the high-energy CP odd observables from the low-energy ones [3] without 
assuming very constraining frameworks for the unmeasurable quantities. Therefore, the 
intent of this work is to clarify the relation between the CP violation accessible to low- 
energy experiments, and the CP violation necessary for leptogenesis, in a phenomenological 
bottom-up perspective, with minimal assumptions about the high scale theory. We just 
assume that the neutrino Yukawa couplings are hierarchical, which is the most natural 
assumption given the observed values of the charged lepton and quark Yukawas. Neutrino 
oscillation data then lead to hierarchical singlet masses. 

In this paper, we aim to answer the phenomenological question of whether the BAU can 
be sensitive to low-energy phases, in the supersymmetric seesaw. We suppose the observed 
BAU is generated via thermal leptogenesis, and enquire whether this restricts the range of 
the phases. A similar issue was investigated by Branco et.al [4], where it was shown that 
for any value of the measurable CP violating phases, a large enough BAU can be produced. 
This statement has been recently confirmed in a study [5] that includes flavour effects [6], 
in the Standard Model seesaw framework. In the present analysis, we want to address 
the question considering flavoured leptogenesis in a supersymmetric scenario, that has the 
interesting feature to potentially add new observables in the lepton sector, through the 
enhancement of flavour and CP violating processes (See eg [7] for a review and references 
on leptonic flavour and CP violation, induced by super symmetry). 

The question we address, and the answer we find, differ from some other analyses 
[8, 9, 10, 11]. As written above, we aim to make few untestable assumptions, and to ask a 
precise phenomenological question: "Is the baryon asymmetry sensitive to PMNS phases?" . 
We find the answer to be no. That is, there is "no correlation" between the BAU and 
PMNS phases, when all the unmeasurables in our scenario are allowed to vary over their 
whole range. To the best of our understanding, Refs. [8, 9, 10, 11] find a correlation 
between the BAU and the PMNS phases because they set unmeasurables (such as phases 
of the "right-handed" neutrinos) to fixed values. 

We define "finding a correlation between Yb and x" to mean "Yg is sensitive to x" . To 
show that the baryon asymmetry Yb is insensitive to (or uncorrelated with) a parameter 
x, we must only show that, for any value of x, we can find a large enough Yb- It would 
be numerically more challenging to show a correlation, because the point distribution in 
scatter plots may reflect the priors on the scanned parameters (see sections 6.4 and 7.2). 
Our definition of correlation differs from that used by [8, 9, 11], and also in [19] (who extract 
correlations from scatter plots). We use our narrow definition because it is parametrisation 
independent. 

Since leptogenesis occurs at a very high-energy scale, a supersymmetric scenario is 
desirable in order to stabilize the hierarchy between the leptogenesis scale and the elec- 
troweak one. However, if supersymmetry exists at all, it must be broken and, in principle, 
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the soft supersymmetry breaking Lagrangian can contain off-diagonal (in flavour space) soft 
terms, that would enhance lepton flavour violating (LFV) processes. These are strongly 
constrained by current experiments; this is the so-called supersymmetric flavour problem. 
In order to avoid it, we focus on the most conservative minimal Supergravity (MSUGRA) 
scenario with real boundary conditions, where the dynamics responsible for supersymmetry 
breaking are flavour blind and all the lepton flavour and CP violation is controlled by the 
neutrino Yukawa couplings. Supersymmetric expectations for LFV [12, 13, 14] and possible 
relations to leptogenesis [7, 11, 15, 16, 17, 18] * and EDMs [19, 20] have been studied by 
many people. 

We perform a scan over the seesaw parameters, looking for those points that give a 
large enough BAU, and where fx — > e^f and one of r — ► ^7 would be seen in upcoming 
experiments. Our analysis is more restrictive than [19], in that we require these branching 
ratios to be "large". The aim is to verify if such experimental inputs imply a preferred 
range of values for the low-energy PMNS phases. We also estimate the contribution to 
the CP violating electron electric dipole moment. A detailed analysis of the MSUGRA 
scenario would require a scan also over the supersymmetric parameters, which is beyond 
the scope of our analysis. 

Due to the large number of unknown parameters, instead of doing a usual grid scan in 
the seesaw parameter space we construct a Markov Chain using a Monte Carlo simulation 
(MCMC — see e.g. [21, 22]). This technique allows to efficiently explore a high-dimension 
parameter space, and we apply it for the first time to the supersymmetric seesaw model 
t. Our work is thus pioneering in the exhaustive scanning of the seesaw parameters, which 
would be otherwise prohibitive without the MCMC technique. 

The paper is organized as follows. In section 2 we introduce the supersymmetric 
seesaw in the MSUGRA scenario and we review the low-energy interactions induced in the 
supersymmetric seesaw model. Section 3 is devoted to thermal leptogenesis with flavour 
effects, and section 4 describes our bottom-up reconstruction procedure. Section 5 gives 
analytic estimates, that complement our numerical analysis, using the MCMC technique, 
which is presented in section 6. We discuss our results in section 7 and conclude in section 
8. 

2. Notation and review 

We consider the superpotential for the leptonic sector in a supersymmetric seesaw model 
[1] with three hierarchical right-handed neutrinos (Mi < M2 < M3): 

M 

W lep = (L L H d )Y e E c + (L L H U )XN C + N c —N c . (2.1) 

In this expression, A, Y e and M are 3x3 matrices, and flavour indices are suppressed. The 
Ll are the supermultiplets containing left-handed lepton fields, E are those containing 

*See ref. [17] for a discussion about when the approximation used in [16] is not valid. 
fSee [23] for a detailed study of the Zee-Babu model of neutrino masses phenomenology using this 
technique. 
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the right-handed charged leptons, while N are the supermultiplets of the right-handed 
singlets. The Majorana mass scale can be taken large 10 9 GeV Mj ^ 10 15 GeV, since 
the corresponding operator is a singlet under the SM gauge group. 

Without loss of generality one can work in the basis where Y e and M are diagonal, so 
that the superpotential gives the following Lagrangian for leptons: 

_ _ 

C = Y ea (t L H*M + (ilK)^*^ + —NW + ... + h.c. (2.2) 

where the parentheses indicate SU(2) contractions and the flavour indices are written ex- 
plicitly. Since supersymmetry is broken, to this Lagrangian we must add the soft SUSY 
breaking terms : 

( BM i 1 

Cssb = E / f / + (—^-i + MVeJl ■ H d e c a + \J% ■ H U N!?) + h.cj (2.3) 

where / collectively represents sfermions. This soft part is written at some high scale Mx 
where, in MSUGRA, the soft masses are universal and the trilinear couplings are propor- 
tional to the corresponding Yukawas. MSUGRA is then characterized by four parameters: 
the scalar (mo) and gaugino (m 1 / 2 ) masses, shared by all of them at the GUT scale; the 
trilinear coupling involving scalars, ao, at the GUT scale; and finally the Higgs vev ratio, 
tan/3. 

In the chosen basis, the neutrino Yukawa matrix is in general not diagonal and complex, 
and can be written as: 

A = vlD x V R (2.4) 

where D\ is diagonal and real. Note that in this basis the neutrino Yukawa matrix is the 
only source of flavour violation in the lepton sector, through the unitary matrices Vl and 
Vr that act respectively on the lepton doublet space and on the right-handed neutrino 
space. These matrices contribute also to CP violation, through six CP violating phases. 
In general, other sources of CP violation appear in the complex neutrino B-term, in the 
scalar mass mo and in the trilinear coupling ao- 

At energies well below the right-handed neutrino mass scale, the effective light neutrino 
majorana mass matrix can be written: 

[m„] = XM~ 1 X T vl = UD U U T . (2.5) 

The first equality shows that the smallness of light neutrino masses is naturally explained 
once the right-handed neutrino mass is set at very high energy, ~ 10 14 GeV (in this 
expression v u = {H u )). In the second equality, D u is a diagonal matrix with real positive 
eigenvalues and U is the PMNS matrix containing the three low-energy CP violating 
phases, the Dirac phase 5 and two Majorana phases a, (3. Those phases are, in general, 
a combination of the 6 phases appearing in the complete theory. We use the standard 
parametrisation: 

e ia C13C12 e %13 S12C13 °- <=~ lh 



U = I e ia (-S12C23 - S23S13C12 e tS ) e lfS (c 23 ci2 - S23S13S12 e tS ) S23C13 j (2-6) 
2 ia (S23S12 - S13C23C12 e tS ) e t/3 (-S23C12 - S13S12C23 e lS ) 
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Present bounds 


Future sensitivity 


BR(r -► /*7) 
BR(r -► 67) 


< 1.2 x 10 -11 

< 6.8 x 10" 8 

< 1.1 x 1(T 7 


10~ 13 (MEG) [251 
10~ 9 (Belle) [26] 


BR(/x -»■ eDeVfj,) 
BR(r -> /xi^i/,-) 
BR(r -> ev e Vp) 


~ 100% 
17.36 ± 0.05% 
17.84 ± 0.05% 





Table 1: Present and predicted bounds on lepton flavour violating processes, and measured branch- 
ing ratios for £ a — > lpv a vp decays. 



Present bounds (e cm) 


Future sensitivity (e cm) 


d e < 1.6 x 10" 27 
dp < 2.8 x 10~ 19 
(-2.2 < d T < 4.5) x 10~ 17 


10~ 29 (Yale group) [27] 
10~ 24 (Muon EDM Collaboration) [28] 



Table 2: Present and anticipated bounds on electric dipole moments. See [7] for a discussion of 
future experiments. 

If we combine the equations (2.4) and (2.5), we can write: 

D v = UWlD x V R D^vSD x VlU*v 2 u = W^D X V R D^V£D X W*V 2 U , (2.7) 

with Vr diagonalizing the inverted right-handed neutrino mass matrix. This relation shows 
that non-zero angles and phases in the unmeasurable right-handed neutrino mixing matrix 
Vr imply non-zero angles and phases in W = VlU, which being in the doublet sector, is 
potentially more accessible. We will use this relation to reconstruct the right-handed sector 
from low energy physics in sec. 4. 

2.1 Low-energy footprints: LFV and EDMs in MSUGRA 

Present bounds on LFV processes, shown in table 1, restrict the size of flavour off-diagonal 
soft terms. This suggests universal soft terms at some high scale Mx, see Eq. (2.3), like 
in the MSUGRA scenario. There are also stringent experimental bounds, as we can see 
in Table (2), on the CP violating electric dipole moments, which point towards very small 
CP phases. To address this "SUSY CP problem" *, we suppose that all the soft breaking 
terms (namely do, mo and right-handed sneutrino B-term), as well as the fi term, are real. 
Even under this extremely conservative assumptions, it is well known that because of RGE 
running from high to low energy scales, the seesaw Yukawa couplings potentially induce 
lepton flavour and CP violating contributions to the soft terms [12, 13, 14]. 

We focus on these neutrino Yukawa coupling contributions to LFV and EDMs, as- 
suming MSUGRA with real boundary conditions at Mx- Additional contributions, arising 
with less restrictive boundary conditions, are unlikely to cancel the ones we discuss, so the 

*See e.g. [24] for an illuminating discussion. 
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upper bounds that will be set if, for instance, no electron EDM is measured by the Yale 
group, will equally apply. Conversely, if an electron EDM is measured above the range that 
we predict, it will prove the existence of a source of CP violation other than the neutrino 
Yukawa phases. 

We are interested in analytic estimates for LFV rates and electric dipole moments. For 
this, we need the flavour-changing and CP violating contributions to the soft masses, that 
arise from the neutrino Yukawa. Following [29], we take the one- loop corrections to the 
flavour off-diagonal doublet slepton masses fnj^ a ^ — > fh^p + Am^ and to the trilinear 
coupling aoA — > aoA(l + Aao) to be: 

Am| Q/3 = ( 3m 2 + a 2 )[C^U ~ (™o + + 2a B)[H} aP , (2.8) 

A(a ) a/3 = -j^ [C^U ~ i^i H U , (2-9) 
for a ^ (3 where the matrices H and are given by: 

H = AA 1 " = VID\V L , (2.10) 

C<"> = A log" » = ^V, log" V^V L . (2.11) 

C^ 1 ) is the leading log contribution, and terms oc iJ arise in the finite part (they could be 
relevant for EDMs). The one loop corrections to the right handed charged slepton mass 
matrix, rh 2 ^^ only contain the charged lepton Yukawa couplings and therefore cannot 
generate off-diagonal entries. These are generated at two loops and, as we will see later, 
they can be relevant for the lepton EDMs. 

At one loop, sparticles generate the dipole operator (where e without subscript is the 
electro-magnetic coupling constant): 

eX^efa^e^F^ + h.c. (2.12) 

which leads to LFV decays (£ a — > £pj), and induces the flavour diagonal anomalous mag- 
netic and electric dipole moments of charged leptons [7]. For a = (3, the anomalous 
magnetic moment is a a = Am ea Re{X aa } and the electric dipole moment is 2lm{X aa }. 

In the mass insertion approximation the observable LFV rates are proportional to 
\fh 2 Laj3 \ 2 oc |C^| 2 and the corresponding branching ratios are of order [13]: 

BR(£ a -> £p 7 ) o? tan 2 /? 2 2 
BR(£ a -+ £pv a vp) ~ G 2 F m% USY '"^ 1 • ) 

a 3 tan 2 /3 (3m 2 , + a 2 ,) 2 2 
~ &F™%VSY ( 47r ) 4 l[CW ' 

where Gjr is the Fermi constant, tan /3 = v u /vd, and msusY is a generic SUSY mass, which 
substitutes for the mixing angles and the function of the loop particle masses. 

An estimate of msusY can be obtained from the data on the anomalous magnetic 
moment of the muon, as suggested in [30]. A 3.3 or 3. 4? deviation from the Standard 
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(116 592 080 ± 63) x 10" 11 in BNK-E821 


5a, = a* xp - a ™ 


(276 ± 81) x 10- 11 [31] 
(275 ± 84) x HT 11 [34] 
(295 ± 88) x lO" 11 [32] 



Table 3: Experimental value and deviation from the SM predictions of the muon anomalous 
magnetic moment. The errors of 5a, are the combination in quadrature of the experimental and 
theoretical ones. 



Model prediction is observed in the anomalous magnetic moment of the muon (in Table 
(3) is given the experimental value of a, and the deviation from the SM prediction [31, 32]). 
We assume it is due to new physics that can also contribute to flavour violation and EDMs. 
In the MSUGRA seesaw scenario that we are considering, the main contribution to a, comes 
from 1-loop diagrams with neutralino or chargino exchange and is given by [33]: 

6a*"" ^ ^ • (2-14) 

8vr sm z weak m z SUSY 

Within this approximation, the observed deviation in the muon anomalous magnetic 
moment only fixes the ratio tan (3/m 2 SUSY ~ 5 10~ 5 GeV~ 2 , so our SUSY masses scale 
with tan (3 as m 2 SUSY = ^(200 GeV) 2 . 

Assuming [30] that the main contribution to the LFV branching ratio is given by 
analogous diagrams involving chargino and neutralino exchange, gives, from equations 
(2.13) and (2.14) with ~ ao ~ msuSY- 

BR{l a ^i pl ) ^ in _ 8 _ . 2 (5a, 



10-»|C^r 773, • (2-15) 



BR(£ a — > £pv a i>p) ' Pl V10 

Since we aim to explore seesaw parameter space, we set the MSUGRA parameters itlq ~ 
a ~ msuSY- 

In our analysis, we aim for values of |C Q| g| 2 that will give fi —> and either of r —> £7 in 
the next round of experiments. We require only one of the r decays, because the other must 
be small to suppress fi — > e^f (recall that we assume the neutrino Yukawas are hierarchical) . 

The neutrino Yukawa corrections to the soft terms can also enhance the predictions of 
the CP violating electric dipole moments. In our discussion we can neglect muon and tau 
EDMs, because the experimental sensitivity on is currently eight orders of magnitude 
weaker than on d e and we expect d,/d e ~ m,/m e . 

There are two potentially important contributions to the charged lepton EDMs induced 
by the neutrino Yukawa couplings. As discussed in [35, 29], the first non-zero contribution 
to the complex, flavour diagonal EDMs arises at two-loop order. The matrices Aao and 
Am 2 in Eq.(2.8) are the available building blocks to make an EDM, which turns out to be 
proportional to the commutator [H, C\. This is the dominant contribution at low tan/?. 
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We follow [29] § to estimate: 

4rv / 9 \ 

4 ~ , . , 5 2 e lm[H C] ee (1.9 10~ n e cm) ~ 10~ 29 — — Im[ff C] e cm , (2.16) 
(4tt) m SUSY \tan/}/ 

where we have used [H,C]/i = 2lm[H C], and the 2/tan/3 arises because we extracted 
m 2 SUSY from the Sa^. 

In the large tan (3 region, it has been shown [36] that a different contribution to the 
EDMs can be the dominant one. This new contribution arises at three loops, and it 
involves the two loop correction to the right handed charged slepton mass matrix Afh 2 E . 
It is proportional to the CP violating quantity: 

D a = Im [((AmlfmAml] aa (2.17) 

where m is the (diagonal) charged lepton mass matrix. Despite being a higher loop order, 
it is typically dominant for tan (5 <^ 10. The two loop expression for Arh 2 E can be found in 
[29]. We approximate this contribution as: 

-e 8a 10m tan/? ImjA^A k m 2 A* m A em ] 

2 (4vr) 7 m 2 SUSY v 2 cos 2 V kh y 1 



where 



, 9N M 2 , M 2 M 2 , 9 M 2 M 2 \ 

F(Ml) = bg log bg — ^ + log 2 log — ^ , (2.19) 



and M x = 3 x 10 16 GeV, M N = M 2 . It gives an electric dipole moment of order: 

o<) /tan/?\ 2 Im[A* fc A Qfc m| A* m A em ] 

a e ~ 10 ecm. 

1 50 / ml 



One comment is in order. Throughout this work, we use the approximated formulae 
(2.15), (2.16) (2.18), where we have set the supersymmetric parameters mo and ao at a 
common msusY scale. Of course these are very rough approximations, but given that 
a detailed analysis of the MSUGRA scenario is beyond the scope of this study, which 
concentrates on the seesaw parameters, it is enough to illustrate our results. 

Notice that, since we normalize the LFV branching ratios to the muon g-2 deviation 
from the SM, there is no enhancement of LFV for large tan (3. The three loop EDM 
contribution (2.18) is enhanced, because it has extra powers of tan/3. 

3. Flavoured thermal leptogenesis 

The observed Baryon Asymmetry of the Universe [37] is: 

= (8.75 ±0.23) x 10~ n (3.1) 



where ubo, n B0 , and so are the number densities of baryons, antibaryons, and entropy, in 
the Universe today. We assume this excess is produced via flavoured thermal leptogenesis [2, 



'[29] finds the same structure as [35, 36], but its result is smaller by one power of a large logarithm. 
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6, 38], through the decays of the lightest singlet neutrino N\ and sneutrino Ni, in the 
thermal plasma at T ~ Mi. The population of N\ and N\ is produced by inverse decays 
and scattering in the plasma. The decays are CP violating and controlled by the neutrino 
Yukawa coupling, thus for hierarchical right-handed (s) neutrinos the CP-asymmetry is 
given by [39]: 



r(iVi -► e a H, £ a h) - r(jVi -► e a H, t a h) 

Y{N 1 £ff, + Y{N 1 -> IE, IK) 
3Mi 



8tto 2 [At A] 



li 



Im{[A] a i[mtA] a i}, (3.2) 



where a specifies the flavour of the (s)lepton doublet in the final state. If the CP violating 
decays are out-of-equilibrium the lepton asymmetry produced can survive and be partially 
converted into a baryon asymmetry through non perturbative SM sphaleron processes [40]. 

In Eq.(3.2) we have intentionally not summed over the flavour index a, because flavours 
can have a role in the evolution of the lepton asymmetry [6]. That is, if a flavour in 
the thermal bath is distinguishable, then the corresponding lepton asymmetry follows an 
independent evolution. This occurs when the charged lepton Yukawa interaction rate 
T( a = 5 x 10 _3 TY^ is faster than the expansion rate H and the singlet inverse decay rate 
Tin ~ e~ m l T YN, where Tjv is the right-handed neutrino decay rate. Since leptogenesis 
takes place at T ~ Mi the mass of the lightest right-handed (s)neutrino tells us if flavour 
effects are important. 

In the MSSM, the charged lepton Yukawas are larger than in the SM: Y a = m Q /(cos f3x 
174GeV), so they come into equilibrium earlier. At very high temperatures T > tan 2 (3 10 12 
GeV ™, the charged lepton yukawa interactions are out of equilibrium {T^ a <C H) and there 
are no flavour effects, so leptogenesis can be studied in one-flavour case. However, as the 
temperature drops, the r interactions come into equilibrium. In the range tan 2 (3 10 9 < 
T < tan 2 [5 10 12 GeV, we have an intermediate two-flavour regime, so that the lepton 
asymmetry produced in the r evolves separately from the lepton asymmetry created in the 
linear combination: 

(3.3) 



ViA + A e ie 



For T < tan 2 j3 10 9 GeV, also the n Yukawa interactions come into chemical equilibrium 
and all the three flavours become distinguishable. 

In all the flavour regimes the baryon to entropy ratio can be written as: 

lOn^v + n^ ^10315C(3) V 
YB -31 s Z^ e ^ ~ 31 4vr 4 Z^ e ^ • V A > 

a 3 a 

The numerical prefactor indicates the fraction of B — L asymmetry converted into a baryon 
asymmetry by sphalerons [41] in the MSSM. The second fraction is the equilibrium den- 
sity of singlet neutrinos and sneutrinos, at T > Mi, divided by the entropy density s. 



'We approximate tan/3 ~ 1/ cos/3 because sin/3 ~ 1 and tan/3 is a more familiar parameter. 
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Numerically, it is of order 4 x 10~ 3 , similar to the non-SUSY case The e aa are the CP 
asymmetries in each flavour (so that a = r, o or a = r, fi, e in the two- or three-flavour 
regimes respectively) and the i] a are the efficiency factors which take into account that 
these CP asymmetries are partially erased by inverse decays and scattering processes. We 
assume the efficiency factors have the same functional form and numerical factors as for 
non-supersymmetric leptogenesis [6]: 



2\A n 



\m r 



-1.16 



+ 



m r 



2m* 



(3.5) 



where we neglect A-matrix [42] factors in our numerical analysis. The rescaled N\ decay 
rate is defined as : 



m 



I A. 



q1 



Mi 



(3.6) 



and in supersymmetry m ^f SSM = m ^ M /y/2 = A-kv^H^/MI ~ 0.78 x 1CT 3 eV **, where H 1 
is the Hubble expansion rate at T = M\. 

Combining equations Eq.(3.4), Eq.(3.2), Eq.(3.5) and Eq.(3.6), we can write the BAU 

as: 

in i ■ ? - i ; 

>t • A] Q }, (3.7) 



10 135Mi ^ - 



31 A^g^v^ 



ml 



where X a = [X] a i/ \/[X^ X]n- Yb is roughly a factor of y/2 larger than in the SM, in the 
limit where m aa > for all flavours. 

Supersymmetric thermal leptogenesis suffers from the so called gravitino problem [43]: 
in a high temperature plasma gravitinos are copiously produced and their late decay can 
jeopardize successful nucleosynthesis (BBN). This gives an upper bound on the reheat 
temperature of the Universe Trh , which constrains the temperature at which leptogenesis 
can take place, and gives an upper bound on the singlet neutrino mass Mi ^ 5Trh [44, 45]. 
However, there is also a lower bound on Mi ^ 10 9 GeV [46] (for hierarchical iVs) to obtain 
a large enough lepton asymmetry. This can be seen from (3.7), where Yb oc Mi. It has 
recently been suggested [47] that this conflict can be avoided by generating the singlet 
masses after reheating. However, we here assume that Mi > 10 9 GeV is fixed before 
reheating. 

There are various ways to obtain Trh ~ 10 9 — 10 10 GeV. If the gravitino is unstable, the 
nucleosynthesis bound leads to very stringent upper bounds on the reheating temperature 
after inflation [48]: T RH & 10 4 - 10 5 GeV for m 3/2 ^ 10 TeV, or T RH & 10 9 - 10 10 GeV for 
777-3/2 > 10 TeV. A sufficiently high reheat temperature is obtained for very heavy gravitinos 
because they decay before BBN. Alternatively, if the gravitino is the stable LSP, a correct 



" The addition of the Ns is compensated by the approximate doubling of the degrees of freedom in the 
plasma : g, = 228.75 for the MSSM. 

"There are factors of 2 for SUSY: defining r D to be the total N decay rate, we have rf, (/5y = 2T S D M . 
So with the definition of eq. (3.6) for in, we have m = 47r^r^ SSM /M 1 2 as opposed to rh = 87r^rf, M /M'{. 

if M /y/2, where m* is the value of in that would give r_o = Hi at T = Mi, and the factor 



SUSY 



So m 

of V2 is because there are approximately twice as many degrees of freedom in the plasma. 
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dark matter relic density can be obtained for T RH ~ 10 9 - 10 10 GeV. In this scenario, one 
must ensure that the decay of the NLSP does not perturb BBN. This can be obtained, 
for instance by choosing the NLSP with care [49] or by having it decay before BBN via 
i?-parity violating interactions [50]. 

We can summarise that a reheat temperature *S 10 9 GeV is difficult but not impossible 
in supersymmetry. So for the purposes of this paper, we will allow Mi < 10 11 GeV. 



4. Reconstructing leptogenesis from low energy observables 

In order to search for a connection between the low-energy observables and leptogenesis, 
we need a parametrisation in which we can input the low energy observables, and then 
compute the BAU. Ideally we want to express the high-energy parameters in terms of 
observables [51]. Therefore, we write the seesaw parameters in terms of operators acting 
on the left-handed space, potentially more accessible: so we chose D v , D\ and Vl (that 
appears in the combination AA^) and Upmns- Within this bottom-up approach, the CP 
violation is now encoded in the three, still unknown, low energy phases of the PMNS matrix 
U, and in the three unknown phases in Vl- We then reconstruct the right-handed neutrino 
parameters in terms of those inputs. 

The matrices D u and Upmns can be determined in low-energy experiments. Through 
neutrino oscillation experiments we can extract the two neutrino mass differences, the 
PMNS matrix mixing angles and, in the future, the Dirac phase [52] (if Nature is kind 
with us). Furthermore, we have an upper bound on light neutrino masses that comes 
from cosmological evaluations [53], Tritium beta decay [54], and neutrinoless double beta 
decay [55]. Observing this last process could prove the Majorana nature of neutrinos and 
put some constraints on the combination of Majorana phases. 

We have seen that in MSUGRA there is an enhancement of lepton flavour violating 
processes due to the neutrino Yukawa couplings. Assuming that these processes can be 
measured in the near future constrains the coefficients [C] a A, see Eq. (2.13), which depend 
on D\ and Vl- We parametrise the Vl matrix as the product of three rotations along the 
three axes, with a phase associated to each rotation: 



I n L „L „L „L e —ip „L „—icr \ 

L 13 L 12 L 13*12 e b 13 e 

ri L C L cAp o — ioj C L f-iL g^icr ri L ri L C L ^ — iu C L C L ^ — ip ^icr ri L C L tD ~ 

L 23*12 e *23 e *13 L 12 e C 23 L 12 6 23 e 6 13 6 12 e e L 13 6 23 e 

L iu> „L ip _ L L L iu _ L iw L _ L L L -ip ia L L 

\ 6 23 e b 12 e 6 13 L 23 L 12 e *23 e L 12 6 13 6 12 L 23 e e C 23 L 13 / 



(4.1) 

^From the bottom-up parameters defined above and using the equation (2.7), we are now 
able to reconstruct the right handed neutrino mass matrix and the Vr matrix appearing 
in the baryon asymmetry: 

M- 1 = VrD^vI = Dl l V L UD u U T VlD- l v- 2 . (4.2) 

In leptogenesis without flavour effects, the BAU is controlled only by the phases of 
Vr, which also contribute to the Upmns i n the parametrisation we use. However, as 
demonstrated in the R matrix parametrisation [56], it is always possible to choose Vl such 



- 11 - 



that the lepton asymmetry e has any value for any value of PMNS phases [4]. So for Yb 
in its observed range, the PMNS phases can be anything, and if we measure values of the 
PMNS phases, Yb can still vanish. In flavoured leptogenesis, the BAU can be written as a 
function of PMNS phases and unmeasurables, but it was shown in [5] that for the Standard 
Model seesaw, Yb is insensitive to the PMNS phases. Relations between low energy CP 
violation and leptogenesis can be obtained by imposing restrictions on the high-scale theory, 
for instance that there are no right-handed phases [8]. 

In the case of MSUGRA, we assume that we will have two more measurable quantities 
in the near future, \i — > e^y and either of r — ► £j. Naively, we do not expect LFV rates to 
add more information on the CP violating phases, because the rates can be used to fix two 
(real) parameters in D\ and Vl- The question is whether the remaining phases and real 
parameters, can always be arranged to generate a large enough BAU. We find the answer 
to be yes. For instance, in the limit of taking only the largest neutrino Yukawa coupling in 
D\, the matrices become proportional to H, and using the parametrisation of the Vl 
matrix given in Eq.(4.1) one can easily see that the CP violating phases of the Vl matrix 
disappear from the LFV branching ratios. 

Besides the LFV processes, the neutrino Yukawa couplings can also contribute to the 
CP violating electric dipole moments. These contributions are expected to be below the 
sensitivity of current experiments [20, 57]. See [57] for a discussion of the impact of EDMs 
on seesaw reconstruction. In our framework with hierarchical Yukawas we expect some 
suppression on this contributions to the EDMs. As we have seen in Section 2.1, for low 
tan/3 the main contribution is proportional to the commutator of the matrices and 
H, see eq. (2.16). Thus in the limit of taking only the largest Yukawa, which implies 
oc H, the commutator is equal to zero. Regarding the large tan/3 regime, although 
the contribution to the EDMs has a different dependence, given in eq. (2.18), it can be 
shown that it also vanishes in this limit. This means that a non-zero contribution will be 
suppressed by mixing angles and a smaller eigenvalue of H. 

5. Analytic Estimates 

If a parametrisation existed, in which one could input the light neutrino mass matrix, the 
neutrino Yukawa couplings that control lepton flavour violation, and the baryon asymmetry, 
then it would be clear that the BAU, and other observables, are all insensitive to each other. 
In this section, we argue that at the minimum values of M\ where leptogenesis works, such 
a parametrisation "approximately" exists. 

We analytically construct a point in parameter space that satisfies our criteria (large 
enough BAU, LFV observable soon), and where the baryon asymmetry is insensitive to 
the PMNS phases. To find the point, we parametrise the seesaw with the parameters 
of the effective Lagrangian relevant to N\ decay. Since the observed light neutrino mass 
matrix is not an input in this parametrisation, one must check that the correct low energy 
observables are obtained. This should occur, in the region of parameter space considered^, 

^This area of parameter space was also found in [58] using a left-handed parametrisation inputting 
W = V L U instead of V L . See also [59]. 
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because the contribution of N\ to the light neutrino mass matrix can be neglected. We 
construct the point for the normal hierarchy and small tan f3; similar constructions are 
possible for the other cases. 

The effective Lagrangian for N\ and Ni, at scale M\ ^ A <C M2, arises from the 
superpotential: 

W Nl = X al L a L H u N^ + ^-iVfiVf + Ka(S (L L a H u )(L^H u ) (5.1) 

where n a p is obtained by integrating out N2 and 7V 3 . It is known [60] that the smallest 
Mi for which leptogenesis (with hierarchical Ni) works, occurs at to* ^ m ^ m so \. So we 
assume that 

-v\ <C m a/ 3, (5.2) 



Mi 

implying that Ni makes negligible contribution to light neutrino observables. We are 
therefore free to tune the A Q is to maximise the baryon asymmetry. 

To obtain a baryon asymmetry Yb — W~ 3 J2 a e aai]a — 8 x 10~ n , we require: 

J2 e <*aV*-8x 10~ 8 • (5.3) 

a 

For tan/3 ~ 2, it is unclear whether the is distinct for leptogenesis purposes. For 
simplicity we assume not, and use two flavours o and r. The efficiency factors rj a are 
maximised to r\ a ~ 1/4 for rh aa = \\ a i\ 2 vl/M! ~ V2m*. Since to ~ 3to*, this is barely in 
the strong washout regime, and (3.5) should be an acceptable approximation. 

We would therefore like to find a point in parameter space, such that Mi ~ 10 9 GeV, 
£oo — £rr — 1-6 x 10~ 7 . Defining X a = X a i/y / ^2 a \X a \ | 2 , equation (3.2) implies that we 
need, for a = o and a = r: 

Im A,M%n^^. (5 . 4) 

I TO3 J Mi 

This means that Ai needs a component along ^3 (the eigenvector of TO3), and, since it 
should also generate mi, it needs a component along u\. It can always be written as: 

Ai = Anni + A 2 i-u 2 + A31U3 , (5.5) 

where {1,2,3} indices indicate the light neutrino mass basis. In the following we take 
A21 = 0, A31 = |A 3 i|e^, I A31 1 > I An I . With equation (2.5), 

I m (A Ql l!^k} = -L- - Im^Xs^ + X 2 31 U a3 )U* a3 } 
1. ma J An r + A3i 1. > 



TO3 J |An| 2 + |A 3 l| 



i.\„i- , + i.v„i- ,J,Q il 



(no sum on a). In the last formula, we drop the terms oc An, which may contain asym- 
metries that cancel in the sum e OQ + e rr . These are not useful to us, because we aim for 
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f]o — Vr — 1/4. For Im {A^j/dAsij 2 + | An | 2 ) k, 1/2, Eq.(5.4) implies that a large enough 
BAU could be produced for M x ~ 3 x 10 9 GeV. 

We now check that we obtain the observed light neutrino mass matrix, even with £, 
the phase of A31, of order tt/4. The light neutrino mass matrix is: 



r 1 A Q lA/3i 2 2 2 

[Ha/3 = M V u + K Q/3 V U = V„ 



\ 2 \ 2 i 

— UiU-l + K 2 U 2 U 2 + (— + K3)U3U 3 ^ (5.6) 



a/3 

where K2 and K3 are the eigenvalues of n. By convention there is no phase on m 3 , so in 
the 2 right-handed neutrino (2RHN) model that generates k, we should put a phase on 
the larger eigenvalue K3. Since X^v^/Mi ~ e t2 ^ x 10~ 3 eV, the phase on k% is very small 
and we neglect it in the following discussion of lepton flavour violation.. It is well known 
[61] that the seesaw mechanism with 2 right-handed neutrinos can reproduce the observed 
light neutrino mass matrix, with mi = 0. In our case, we assume that A 2 and A3 give the 
observed m 2 , and m 3 up to (negligeable) corrections due to Aq of order 10 -3 eV. mi arises 
due to N±. 

In the 2RHN model, there is less freedom to tune the LFV branching ratios [62] 
than in the seesaw with three iVj. So as a last step, we check that we can obtain LFV 
branching ratios just below the current sensitivity. The 2RHN model can be conveniently 
parametrised with D K , the 3x2 Upmns matrix, the 2x2 unitary matrix W = VlU, and 
the eigenvalues A 2 and A 3 of A (matrices in the 2RHN subspace are denoted by hats) . A is 
a 2 x 2 sub-matrix of A, obtained by expressing the 3x3 Yukawa matrix in the eigenbases 
of the heavy and light neutrinos, and dropping the first row and column, corresponding to 
v\ and N\. It is straightforward to verify that [Vl]3 C ~ 10 -3 can be obtained by taking 
tan^u/ ~ 513/(013512), where 6\v is the rotation angle in W and 6^ are from Upmns- 
Choosing A 2 , the smaller eigenvalue of A, to be ~ .06, ensures that BR([x — > ej) is small 
enough. We can simultaneously take A3 ~ 1 and obtain \Vl\zt ~ [Vl]3/* ~ 1, which allows 
BR(t -> ni) ~ 10" 8 . The resulting masses of N 2 , N 3 are ~ 10 12 , 10 15 GeV. 

Our MCMC has some difficulties in finding the analytic points. We imagine this to 
be because they are "fine-tuned" in the parametrisation used by the MCMC. The amount 
of tuning required in the angles of Vl, to obtain the desired {Aji}, can be estimated by 
taking logarithmic derivatives. In Appendix A, we find a fine-tuning of order: 

™ 2 01 (5.7) 



mp 13 

where 9ij are the Upmns phases, and we optimistically assumed #13 ~ .1. These points at 
Mi ^ 10 10 GeV with m <; 10~ 3 eV, were also not found in the analysis of [19]. 



6. MCMC 

In this section we describe our numerical analysis. In order to verify if the baryon asymme- 
try of the universe is sensitive to the low energy PMNS phases, we perform a scan over the 
neutrino sector parameters aiming for those points compatible with the measured baryon 
asymmetry and the bound on the reheating temperature, that have large enough LFV 
branching ratios to be seen in the next experiments. 
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Light neutrino best fit values 



Am 2 sol = (7.60 ± 0.20) x 1(T 5 eV 2 
Am 2 atm \ = (2.40 ± 0.15) x 10" 3 eV 2 
sin 6 2 sol = 0.320 ±0.023 
sin 2 ,„, = 0.500 ±0.063 



Table 4: The best fit values of the light neutrino parameters and their la errors [63]. 

Using the bottom-up parametrisation of the seesaw defined by the Vl, D\, D u and U 
matrices, our parameter space consists of the 14 variables displayed in Table 5. We take as 
an experimental input the best fit values of the light neutrino mass differences and of the 
solar and atmospheric mixing angles, Table 4. With respect to the SUSY parameters, we 
choose two different regimes for tan /3, equal to 2 or 50, while the msusY scale is deduced 
from the data on the anomalous magnetic moment, see section 2.1. 

Due to the large number of parameters it would prohibitive to consider a usual grid 
scan. Thus, we choose to explore our parameter space by a Markov Chain Monte Carlo 
that behaves much more efficiently, and has been already successfully employed in other 
analyses [64]. 

6.1 Bayesian inference 

Given a model with free parameters X = {xi, . . . ,x n } and a set of derived parameters 
£(X), for an experimental data set d, the central quantity to be estimated is the posterior 
distribution P(X\d), which defines the probability associated to a specific model, given the 
data set d. Following the Bayes theorem, it can be written as: 



where C(d\X) is the well known likelihood, that is the probability of reproducing the 
data set d from a given model X, ir(X) is the prior density function, which encodes our 
knowledge about the model, and P(d) = J C(d\^(X))n(X)dX is an overall normalization 
neglected in the following. In the case of flat priors: 



the posterior distribution reduces to the likelihood distribution in the allowed parameter 
space. 

The main feature of the Markov chains is that they are able to reproduce a specific 
target distribution we are interested in, in our case the posterior distribution, through 
a fast random walk over the parameter space. The Markov chain is an ordered se- 
quence of points Xi with a transition probability W(Xi + i\Xi) from the i — th point to 
the next one. The first point Xq is randomly chosen with prior probability tt(X). Then 
a new point is proposed by a proposal distribution Q(Xi + \\Xi) and accepted with prob- 
ability A(Xi + i\Xi). The transition probability assigned to each point is then given by 



P(X\d) 



c(d\t(x)Mx) 
p(d) 



(6.1) 




(6.2) 
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W{Xi + i\Xi) = Q(Xi + i\Xi)A(Xi + i\Xi). Given a target distribution P(X), if the following 
detailed balance condition: 

W(X k \Xj)P(Xj) = W(Xj\X k )P(X k ) (6.3) 

is satisfied for any j, k, then the points Xi are distributed according to the target distribu- 
tion. For a more detailed discussion see [21, 22]. 

6.2 The Metropolis-Hastings algorithm 

In order to generate the MCMC with a final posterior distribution (6.1), we use the 
Metropolis-Hastings algorithm. In the following, we briefly recall how the algorithm be- 
haves, but the discussion is done in terms of the likelihood, instead of the posterior distri- 
bution, since we assume flat priors on our parameter space, see eq. (6.2). 

Let X be the parameter set we want to scan, and C(X) our likelihood function, the 
target distribution. /,From a given point in the chain with likelihood £(JQ), a new point 
X new with likelihood C(X new ) is randomly selected by a gaussian proposal distribution 
Q(X new , Xi) centered in Xi and having width e. This last quantity e controls the step size 
of the random walk. The new point is surely added to the chain if it has a bigger likelihood, 
otherwise the chain adds the new point with probability C(X new ) / C(Xi) . So the value of 
the next point Xi + \ in the chain is determined by: 

x ^_\x new with probability mm[A(X new ,Xi),l] 
) Xi with probability 1 — mm[A(X new , Xi), 1] 

where A(X new , Xi) is the acceptance probability: 

A(X new ,X l ) = ^^. (6.5) 

Given this acceptance distribution and using the symmetry of our proposal distribution 
Q(Xi,Xi) under the exchange I <-> i, it is straightforward to see that the detailed balance 
condition 6.3 is satisfied for the likelihood C(X) as target distribution. This implies that 
when the chain has reached the equilibrium, after a sufficiently long run, our sample is 
independent of the initial point and distributed according to C(X). 

In order to arrive at the equilibrium in a reasonable amount of time, the step scale e 
of our random walk must be accurately chosen. Indeed, if we define the acceptance rate as 
the number of points accepted over the number of points proposed, a too big step e implies 
a too low acceptance rate, so that our Markov Chain never advances, while a too small e 
and, so, a very large acceptance ratio, implies that our chain needs a very large time to 
scan all the space. It has been suggested that e must be chosen according to an optimal 
acceptance rate between 20% and 50%. However, in order to ensure the detailed balance 
condition, e cannot change during the run of the chain, thus, it is set by our program in a 
burn- in period. 

A valid statistical inference from the numerical sample relies on the assumption that 
the points are distributed according to the target distribution. The first points of the 
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Free parameters 


Allowed range [X min , A" ma:E ] 




A 2 /Ai ~ A3/A2 ~ 30 


A2/A1 ~ 100, A3/A2 ~ 50 


login A"? 


[-0.3,0.3] 


[-0.5,0.5] 


login A? 


[-1.77,-1.17] 


[-2.2,-1.2] 


login Ai 


[-3.25,-2.65] 


[-4.2,-3.2] 


iU oiov /n/ l/ c v ^ 


[-6,-3] 


logio^ 


[-4, log 10 it] 


p, CJ, 0" 


[0,tt] 


#13 


[0..0.2] 




[0,7T] 


a, /3 


[0,tt/2] 



Table 5: Allowed parameter space, so that the uniform prior on each parameter is defined as in 
eq.(6.2). 



chain are arbitrarily chosen and the chain needs a burn-in period to converge to the target 
distribution. The length of the burn-in strongly depends on the intrinsic properties of the 
chain and cannot be set a priori. It changes according to the complexity of the model, 
to the target distribution, and the efficiency of the proposal distribution employed. Once 
the chain has reached the equilibrium the first burn-in points must be discarded to ensure 
the independence of the chain from the initial conditions. Nevertheless, as we will see in 
section 6.4, even following the procedure above, it can be a delicate issue to determine if a 
chain has really converged. 

6.3 The seesaw sample 

In our work the free variables X are given by the 14 seesaw parameters, with uniform 
priors, Eq. 6.2, on the allowed range of parameter space (see Table 5). The choice of a 
logarithmic scale on some unknown parameters allows us to scan with the same probability 
different orders of magnitude. We analyze models with two different hierarchies in the 
neutrino Yukawas, so that, for a A3 ~ 1 we impose A2/A1 ~ A3/A2 ~ 30 or A2/A1 ~ 100 
and A3/A2 ~ 50. The lightest neutrino mass is allowed to vary between three orders of 
magnitude 10~ 6 < m\ < 10~ 3 eV and the #13 mixing angle within its 3d range, < #13 < 
0.2 rad. The Vl mixing angles can vary over 4 orders of magnitude, with maximum value 
it. All the CP violating phases, those of the Vl matrix indicated by p, uj and a and the 
Dirac and Majorana phases 5, a and (3, are allowed to vary on all their definition range: 
[0, it/2] for the Majorana phases and [0, it] for the others (this avoids degeneracies). 

The idea is, now, to generate a sample of points in our parameter space that provide 
enough BAU, give LFV rates big enough to be seen in the next generation of experiments, 
and also have an M\ light enough to avoid the gravitino problem. We then define our set 
of derived parameters £(X) as in Table 6 and we associate to them a multivariate gaussian 
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likelihood with uncorrelated errors: 



£{£exp\£,) — 



1 



(6.6) 



(2vr) 1 /2^m/2 



Where m = 4 is the dimension of the derived parameter set. The centre values ^ exp are the 
best fit values and 1Z is an m x m error matrix, in this case diagonal, since we assume no 
correlation between the errors. As we can see in Table 6, the BAU is set to its experimental 
value, while the LFV rates are set to be one order of magnitude below the present bounds, 
and the expected value of lightest heavy neutrino mass Mi ~ 10 9 GeV is set to escape 
the gravitino problem. The branching ratio of LFV r decays is given in terms of the 
combination BR(t — > ej) + BR(t — > fj,j) = BR ra , since one of them is suppressed to 
respect the stringent bound from BR(/j, — ► ej) (we assume hierarchical yukawas). 

For each point of the chain, the lepton flavour violating branching ratios are es- 
timated with equation Eq.(2.15), while Yb is computed after the reconstruction of the 
right neutrino mass, see Eq.(4.2), using Eq.(3.4) in the flavour regime is in act at the 
temperatures we consider. We recall that the temperature at which leptogenesis takes 
place is of the same order of the reconstructed right-handed neutrino mass. Depending 
on the value of tan f3, the range of temperatures at which the flavour regimes have a role 
changes. As we already mentioned in Section 3: for small tan (3, in the temperature range 
10 9 GeV < T < 10 12 GeV the r flavour is in equilibrium and the two flavour regime is 
in order; while for T < 10 9 GeV \i are also in equilibrium and the three flavours are dis- 
tinguishable. Since we aim for values of M\ ~ 10 9 GeV if we consider a small value of 
tan (5 our program takes into account that the BAU can be produced in both two or three 
flavour regimes. For very large tan/3, instead, already for T < 10 12 GeV r and fj, are in 
equilibrium, thus the three flavour regime always takes place. 

In the case of steeper yukawa hierarchy, in agreement with our analytical estimate, we 
enlarge our set of derived parameters and maximise the rescaled N\ decay rate to rh ~ 1CT 3 
eV and the heaviest right-handed neutrino masses to M2 ~ 10 12 GeV and M3 ~ 3 10 14 



All the points that do not respect the present bounds on LFV, do not have large enough 
baryon asymmetry or have M\ > 10 11 GeV, have a null likelihood. We assume that the 
largest uncertainty on the baryon asymmetry comes from our calculation, so we allow Yg 
to be as small as 4 10~ n . Those points having one of the RH neutrino masses above the 
Mgut ~ 3 10 16 GeV scale have a null likelihood too, since in that case the equations we 
use for the evaluation of LFV processes do not apply. 

6.4 Convergence 

Convergence of the chain ensures the sample is distributed according to the target dis- 
tribution and thus allows to be confident of its statistical information. The question we 
want to answer in this paper, however, does not require a statistical interpretation of the 
sample. Here we only aim to show that, for any value of the low energy phases, the unmea- 
surable high energy parameters can be rearranged to obtain the right baryon asymmetry. 
Therefore a careful diagnostic of the convergence is not a priority. Nevertheless, we briefly 



GeV. 
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Derived parameters £(X) 


Cexp ± 0" 


Y B 


(8.75 ±0.23) 10" 11 


logio BR(ji- > &y) 


-13 ±0.1 


log 10 BR(T->l 7 ) 


-9 ±0.1 


log 10 (M 1 /GeV) 


-9 ±0.1 



Table 6: Best values and errors for the derived parameters £(X) we want to maximize. 

discuss it in this section since it is an important issue that can help the reader to have 
a better overview on our results. Our sample, indeed, has some typical features that can 
make difficult to check if the chain has reached the target distribution. 

As a rudimentary attempt, in our analysis we use the simplest and straightforward 
approach. We run different chains starting from different values and compare the behaviour 
of the free parameters, once the chains have converged they should move around the same 
limiting values. However, this method can be inadequate in case of poor mixing, i.e. when 
the chains are trapped in a region of low probability relative to the maximum of the target 
distribution. This happens in models with strongly correlated variables, when the proposal 
distribution does not efficiently escape this region. Therefore, it can be an issue for our 
numerical analysis, when, as mentioned in section 5, we look for a fine-tuned region with a 
large baryon asymmetry and low M\. We can understand the poor mixing situation if we 
imagine a landscape on the parameter space corresponding to the target distribution, with 
some broad hills and a tall but very thin peak at the maximum of the target distribution. 
In that case, the step of the chain can be optimized to efficiently scan all the space but, if 
its size is larger than the width of the peak, it can easily miss it. 

In case of strongly correlated variables it can also happen that the region to be scanned 
is mainly a plane, that is with almost null likelihoods. This is the case of our sample, where 
we expect a large region with null or almost null likelihood, for all those points that do 
not have large enough baryon asymmetry, low M\ or do not respect the bounds on LFV. 
In this context, if a gaussian-like proposal distribution, as in our sample, is employed, the 
choice of the starting point becomes important to allow the chain to advance. Indeed, if the 
initial value is surrounded by points with null likelihood (and so null acceptance rate) and 
its distance from the interesting region is much larger than the step of the random walk, 
the chain cannot move from this point, since it always finds points with null likelihood. 
On the other side, if the chain starts in a region which is a reasonable fit to the data, it 
advances. Discarding the first points of the chain can ensure independence of the chain of 
the initial conditions inside the interesting region however, if this region is well separated 
from another interesting region, the chain has almost null probability to find the second 
one. 

In order to perform a valid statistical analysis, more sophisticated methods should 
be employed to decide if the chain has converged. In literature many studies exist on 
convergence criterion that help to check the mixing of the sample and are based on the 
similarity of the resulting sampling densities of input parameters from different chains. An 
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example can be found in [65] and [66]. 
6.5 Run details 

In this subsection we explain the details of our MCMC run. The parameter space we 
scan is very large if compared to the derived variables and, in addition, we expect a strong 
correlation between the evaluated baryon asymmetry and the lightest right-handed neutrino 
mass, see eq. 3.7. Thus, since we expect a sample with poor mixing, as discussed in section 
6.4, we first look for an initial point which is a reasonable fit to our observables. This 
procedure is done running previous shorter chains without imposing null likelihoods to the 
not interesting points. Once a wide enough set of interesting starting points is found, we 
start running the chains. 

All the simulations we present are performed by running 5 chains with 10 6 points each. 
As explained before, during the first burn-in iterations, the scale of the random walk e is 
varied until the acceptance rate of points is between the optimal range 20% and 50%. This 
usually takes much less than 3 10 3 iterations. When the optimal acceptance rate is reached, 
the scale e is fixed during the rest of the run. The chains are then added together after 
having discarded the first 10 5 points, corresponding to the burn-in period, in order to give 
enough time to the chain to converge. As discussed above, this procedure should eliminate 
the dependence on the initial point inside the interesting region, but is only a first attempt 
to ensure the sample has reached equilibrium. We run simulations for both normal and 
inverted hierarchy, in the two cases of small and large tan (3. 

7. Discussion 
7.1 Assumptions 

We assume a three generation type I seesaw with a hierarchical neutrino Yukawa matrix. 
We require that this model produces the baryon asymmetry via flavoured thermal lepto- 
genesis, and induces the observed light neutrino mass matrix. This model has a hierarchy 
problem, so we include supersymmetry. 

We make a number of approximations and assumptions in supersymmetrising the see- 
saw. First, we use real and universal soft terms at some high scale, above the masses of 
the singlet neutrinos. In this restrictive model, the only contributions to flavour off-diagonal 
elements of the slepton mass 2 matrix = [rfi 2 ] Q/ 3, arise due to Renormalisation Group run- 
ning. Second, we use simple leading log estimates for the off-diagonals [m 2 ] a/ 3. Third, we 
estimate the SUSY contributions to the dimension five dipole operator (see Eq.2.12) using 
simple formulae of dimensional analysis (see equations (2. 15), (2. 16), (2.18)). This operator 
induces flavour diagonal electric and magnetic dipole moments, and the flavour changing 
decays i a — > £37. We assume the (g — 2) M anomaly is due to supersymmetry, and use it 
to "normalise" the dipole operator. This implies that our SUSY masses scale with tan/3: 
m 'susY = -^y^(200 GeV) 2 . We imagine that there is an uncertainty ~ 10 in our estimates 
of electric dipole moments and i a — > ^37 decays rates, due to mixing angles and sparticle 
mass differences. 
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Our first approximation, of universai soft terms, seems contrary to our phenomenologi- 
cal perspective: the RG-induced contributions to [m 2 ] a/ 3 can be interpreted as lower bounds 
on the mass 2 matrix elements. However, we neglect other contributions, and require that 
the RG induced flavour-violating mass terms are oc Ce$ (see eq. (2.11)), give detectable 
rates for ji — ► ej and r — > £7 in upcoming experiments. Realistically, measuring \i — ► 
mediated by sleptons might allow to determine rr? eji , but does not determine the seesaw 
model parameters C^} . This model dependence is compatible with our phenomenological 
approach, because our result is negative: we say that even if we could determine ci]} , the 
baryon asymmetry is insensitive to the PMNS phases. 

In our numerical analysis we sample the lightest neutrino mass mi and the PMNS 
mixing angle #13, but these two low energy parameters could be eventually measured. In 
this case our simulations should be reconsidered. However, from the analytical estimates, 
we do not expect that fixing these parameters will change our conclusions. 

7.2 Method 

We explore the seesaw parameter space with a Monte Carlo Markov Chain, for two reasons. 
First, an MCMC is more efficient than a grid scan for multi-dimensional parameter space. 
It is essentially a programme for exploring hilltops in the dark. Since the programme likes 
to step up and is reluctant to step down, it takes most of its steps in the most probable 
areas of parameter space. 

The second potential advantage of a MCMC, is that it could make the results less de- 
pendent on the priors, that is, the choice of seesaw parametrisation, and of the distribution 
of points. The results of parameter space scans are often presented as scatter plots, and 
it is difficult to not interpret the point distribution as probability. However, the density of 
points in the scatter plots depends not only on what the model predicts, but also on the 
distribution of input points. For this reason, seesaw scans using different parametrisations 
can distribute points differently in scatter plots. For example, if a model parameter such as 
a Yukawa can vary between and 1, the results will be different depending on whether the 
Yukawa is 0(1) (take points uniformly distributed between and 1) or can vary by orders 
of magnitude (take the exponential of a variable uniformly distributed between — n and 0). 
We had hoped that an MCMC could improve this, because a converged MCMC distributes 
points in parameter space according to a likelihood function. However, in practise there 
are various difficulties. 

The prior on the seesaw model parameter space matters, because the MCMC takes 
steps of some size in each parameter: broad hilltops are easier to find than sharp peaks. 
As discussed in [66], this can be addressed by describing the model with parameters that 
match closely to physical observables. For this reason we parametrise the seesaw in terms 
of the diagonal singlet mass matrix Dm, the light neutrino mass matrix m = UD U U T ', and 
the neutrino Yukawa matrix AA+ = vId\Vl- These are related to low energy observables, 
because AA^ controls the RG contributions to the slepton mass matrix. We take the priors 
for our inputs as given in Table 5. However, the baryon asymmetry and the mass M\ 
belong to the "right-handed" sector, so are complicated functions of the "left-handed" 
input parameters. The bridge between the LH and RH sector is the Yukawa matrix, whose 
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Figure I: Density of "successful" points, as a function of the lightest right-handed neutrino mass 
Mi and rescaled decay rate m, assuming A3 ~ 1 and A2/A1 ~ A3/A2 ~ 30, for two different 
simulations: NH and tan/3 = 50 (left), and IH and tan/3 = 2 (right). "Successful" points have 
Yb > 4 10~ n , and BR(fi — » ej) and BR(t — > £j) an order of magnitude below the current bounds. 
See section 6.3. 



hierarchies may strongly distort the MCMC step size. To obtain a large enough baryon 
asymmetry for Mi ~ 10 9 GeV requires careful tuning in the "right-handed" space, and our 
MCMC has difficulty to find these points. This is related to a second, practical problem, 
that there are many more parameters than observables, so the space to explore is big, but 
the peaks with enough baryon asymmetry and small enough Mi are rare. It is difficult to 
ensure that the MCMC has found all the peaks, as is discussed in section 6.4. 

In section 5, we find analytically an area of parameter space that satisfies our con- 
straints, but where the baryon asymmetry is insensitive to PMNS phases. This area 
corresponds to the limit where N\ makes a negligible contribution to the light neutrino 
mass matrix. In this area, the seesaw model can be conveniently parametrised with the 
interactions of the effective theory at Mi, and it is straightforward to tune the coupling 
constants to fit the light neutrino mass matrix, LFV rates, and the baryon asymmetry. 

7.3 Results 

The aim of our analysis was to verify if a preferred range of values for PMNS phases S, a 
and (3 can be predicted, once low energy neutrino oscillation data, a large enough BAU, 
and LFV processes within the sensitivity of future experiments are requirements of the 
model. 

In Fig. I, we show the distribution, as a function of the singlet mass Mi and the 
(rescaled) decay rate mi, of the successful points for a yukawa hierarchy A2/A1 ~ A3/A2 ~ 
30, with A 3 ~ 1. 

With the parametrisation described in section 6.3, the MCMC easily finds larger values 
of Mi and m, than the "tuned" points found analytically in Section 5. This preference for 
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Figure II: Density of "successful" points, as a function of the baryon asymmetry and the lightest 
right-handed neutrino mass, assuming A3 ~ 1 and A2/A1 ~ A3/A2 ~ 30, for two different simula- 
tions: NH and tan/3 = 50 (left), and IH and tan/3 = 2 (right). "Successful" points are defined as 
for Figure I. 

larger M\ is expected, because the baryon asymmetry and right-handed neutrino masses 
are strongly correlated, see Fig. II and eqn (3.2). 

Nonetheless, as illustrated in Fig. Ill, the MCMC succeeded in finding points at lower 
Mi, with a steeper M hierarchy in the yukawas A3 ~ 1, A2/A1 ~ 100 and A3/A2 ~ 50. The 
difficulties of finding these tuned points are discussed in section 6.4. 

The importance of the ~ 2 decrease in M\ and to, at the tuned points, is unclear to us: 
the cosmological bound is on Trh, rather than M\. Since in strong washout, an equilibrium 
population of N\ can be generated for Trh *5 Mt/5, the points found by the MCMC at 
M\ ~ 10 10 GeV, could perhaps generate the BAU at the same Trh as the analytic points. 
In any case, we see in Fig. II that the fraction of points with big enough Yb is very sensitive 
to Mi, and therefore to details of the complicated reheating/preheating process. 

In Fig. IV, we show density plots of the points resulting from our Markov Chains, 
corresponding to the the yukawa hierarchy A2/A1 ~ A3/A2 ~ 30, with A3 ~ 1, for normal 
hierarchy (NH) of the light neutrino masses and tan/3 = 2, and for inverse hierarchy (IH) 
and tan/3 = 50. In Fig. VII (plot on the left) we show a density plot in the 5 — (3 
plane for tan/3 = 2, NH and the steeper hierarchy A2/A1 ~ 100, A3/A2 ~ 50 and A3 ~ 1. 
From those plots we see that, for any value of the phases 5, a and j3 our conditions are 
satisfied. The analytic results of Section 5 agree with this. Thus, we can conclude that the 
baryon asymmetry of the universe is insensitive to the low energy PMNS phases, even in 
the "best case" where we see MSUGRA-mediated lepton flavour violating processes. For 
completeness we also show correlation plots between the generated BAU and the three low 
energy phases in Fig.V. The low energy observables do not depend on tan (3, because we 
assume the (g — 2)^ discrepancy is due to slepton loops, and we use it to normalise the LFV 

^The smallest yukawa must be small enough to ensure rh ~ m„. 
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Figure III: Density of "successful" points, as a function of the lightest right-handed neutrino mass 
Mi and rescaled decay rate rh, on the left-side, and between the baryon asymmetry and the lightest 
right-handed neutrino mass, on the right-side. We assume here A3 ~ 1 and A2/A1 ~ 100 and 
A3/A2 ~ 50, for a NH in the light neutrinos and tan/3 = 2. "Successful" points are defined as for 
Figure I. 



rates (see Eqn. 2.15). On the contrary, the value of tan (3 is relevant in leptogenesis because 
it changes the number of distinguishable flavours. However, as we can see comparing plots 
for small/large tan/3, the value of tan/3 does not change our conclusions. 

In Figs. VI and VII (plot on the right), we plot the contribution to the electric dipole 
moment of the electron, arising in the MSUGRA seesaw with real soft parameters at the 
high scale. For both low and large tan /3, points from our MCMC generate an electron 
EDM £ 10- 30 ecm. This agrees with the results of [29, 35, 20]. 

8. Summary 

The aim of this work was to study whether the baryon asymmetry produced by thermal 
leptogenesis was sensitive to the "low energy" phases present in the leptonic mixing matrix 
Upmns- We considered the three generation type-I supersymmetric seesaw model, in 
the framework of MSUGRA with real soft parameters at the GUT scale, and required 
that it reproduces low energy neutrino oscillation data, generates a large enough baryon 
asymmetry of the Universe via flavoured leptogenesis and induces lepton flavour violating 
rates within a few orders of magnitude of current bounds. We then enquired whether a 
preferred range for the low energy PMNS phases 5 and (5 can be predicted. 

We used a "left-handed" bottom-up parametrisation of the seesaw. Our parameter 
space scan was performed by a Monte Carlo Markov Chain (MCMC), which allows to 
efficiently explore high-dimensional spaces. It prefers to find the right-handed neutrino 
mass Mi 10 10 GeV, but can also find successful points with a smaller M\ if it takes small 
steps in the relevant area of parameter space. In this area, we can also show analytically 
that the baryon asymmetry is insensitive to the PMNS phases. 
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Figure IV: Density plots in the plane of the low energy phases S — a and 5 — f3 in models with 
A3 ~ 1 and A2/A1 ~ A3/A2 ~ 30. Upper plots correspond to a simulation with NH and tan/3 = 50, 
and lower plots to IH and tan/3 = 2. "Successful" points are defined as for Figure I. 



We have checked that there is no correlation between successful leptogenesis and the 
low energy CP phases. That is: for any value of the low energy phases, the unmeasurable 
high energy parameters and the still unmeasured mi and #13 can be arranged in order to 
have successful leptogenesis and LFV rates in the next round of experiments. The analytic 
estimates indicate that this result will still be true even if mi and #13 are measured and 
fixed to their experimental values. Finally, we have estimated, for each point in our chains, 
the contribution of the complex neutrino Yukawa couplings to the electric dipole moment 
of the electron. As expected, we find it to be ^ 10 -30 ecm, just beyond the reach of next 
generation experiments. 
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Figure V: Density of "successful" points, as a function of the BAU and the low energy phases in 
models with A3 ~ 1 and A2/A1 ~ A3/A2 ~ 30. Upper plots correspond to a simulation with NH and 
tan/3 = 50, and the lower plots to IH and tan/3 = 2. "Successful" points are defined as for Figure 
I. 
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Figure VI: Density of "successful" points, as a function of the baryon asymmetry and the electron 
EDM generated by neutrino yukawas in models with A3 ~ 1 and A2/A1 ~ A3/A2 ~ 30. The left 
panel corresponds to a simulation with NH and tan /3 = 50, and the right panel to IH and tan (3 = 2. 
"Successful" points are defined as for Figure I. 
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Figure VII: Density plot in S — (3 plane and correlation between the baryon asymmetry and 
the electron EDM generated by neutrino yukawas. We assume here A3 ~ 1, A2/A1 ~ 100 and 
A3/A2 ~ 50, for a NH in the light neutrinos and tan/3 = 2. "Successful" points are defined as for 
Figure I. 



A. Fine tuning of the analytic points 

In this Appendix, we estimate the fine-tuning of the points discussed in section 5, with 
respect to the parametrisation of section 4, which is used by the MCMC. 

We do this in two steps. First, in the parametrisation of section 5, we estimate the 3x3 
matrix which diagonalises m in the basis where AA^ is diagonal. Approximating 

this diagonal Yukawa basis to be the one where AA^ is diagonal, we obtain: 



1 









w 



(A.l) 



where [£W]t is the small rotations that rediagonalise m = (Ay + D K )v^, and Ay = 
AiiAji/Mi. If is parametrised as in eqn (4.1) (but neglecting phases for simplicity), we 
find 
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(A.4) 



To obtain A21 negligeable compared to A31 in eqn (5.5), requires no particular tuning of 
0J2 an d @i3 with respect to A21 andAsi. 

The second step is to estimate the tuning required to obtain small angles 9^ and 6^ 
111 With parametrised as in eqn (4.1), this happens if the angles of Vl 
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satisfy djj ~ Oij (for i,j = 12, 13). So the "tuning" required in 6*f 2 and 6± 3 to obtain small 



-^--^ - — o (A.5) 

^2^3 ™3#13 1 ' 

This implies that #f 3 must be tuned against #13 to obtain 6^ ~ .01. If instead 6*13 ^ .01, 
there is no particular tuning of 6^, and the tuning of 6^ with respect to d\ 2 is or order 
ra/771,3. 
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